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Abstract 

Critical and compensation properties of a mixed spin-1 and spin-3/2 Ising fer- 
rimagnet on a square lattice are investigated by standard and histogram Monte 
Carlo simulations. The critical temperature is studied as a function of a single- 
ion anisotropy strength. The second order of the phase transition is established by 
finite-size scaling for the entire boundary. Some previously obtained results, such as 
a tricritical point, predicted by the mean field theory (MFT) and effective field the- 
ory (EFT), or a first-order transition line separating two different ordered phases, 
obtained by the cluster variational theory (CVT), are deemed artifacts of the respec- 
tive approximations. So is a reentrant phenomenon produced by CVT. Nevertheless, 
the multicompensation behavior predicted by MFT and EFT was confirmed. 
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1 Introduction 



There have been a number of theoretical studies of mixed-spin Ising fer- 
rimagnets as simple models of certain types of molecular-based magnetic ma- 
terials [IHS]. Besides other interesting properties, such as appearance of a 
multicritical behavior, they attract attention due to so-called compensation 
behavior with possible technological applications. Such a phenomenon occurs 
at a compensation temperature, i.e., the temperature below the critical point 
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at which the sublattice magnetizations completely cancel out and the total 
magnetization changes sign. Generally speaking, from the previous studies it 
can be concluded that the higher values of spins and the more complex lattice 
topologies used the richer behavior can be expected. However, the increas- 
ing complexity usually requires increasing simplifications in the approaches 
for the problem to be manageable. Thus nonperturbative approaches, such as 
exact [GHTO] or Monte Carlo (MC) [TTHTi] . are so far limited to the simplest 
cases of either the smallest spin values (i.e., mixed spin-1/2 and spin-1) or the 
lattice topology (e.g., honeycomb or Bethe lattices). As a consequence, the 
behavior of the mixed spin-1/2 and spin-1 Ising system is rather well under- 
stood, while there are still disagreements among different theoretical investi- 
gations of the mixed-spin systems with higher spin values, including the spin-1 
and spin-3/2 case. The disagreements arise from the fact that due to higher 
complexity mostly different approximative schemes with various degrees of 
approximation have been employed, such as mean field theory (MFT) [T5] . 
effective field theory (EFT) [TGHIH] and cluster variational theory (CVT) |19j . 
These approximative approaches have been previously shown to produce some 
artifacts, such as a tricritical point, predicted by MFT |20] and EFT |2T] for 
the mixed spin-1/2 and spin-1 system on a square lattice, that were not repro- 
duced either in numerical transfer matrix |12J or MC studies [II1IT21II1]. To 
our knowledge, the results of these approximative studies for the mixed spin-1 
and spin-3/2 system have so far been verified by MC simulations only for the 
case of a simple cubic lattice [22]. However, for this particular case, there were 
basically no qualitative disagreements among the conclusions drawn from the 
respective studies and the MC results only confirmed the previously obtained 
results. On the other hand, for the cases of honeycomb and square lattices, 
there are qualitative differences in the results, which have not been resolved 
yet. 

The objective of this study is to focus on the case of the mixed spin-1 and 
spin-3/2 Ising system with a uniform single-ion anisotropy on a square lattice, 
for which the differences between the MFT, EFT and CVT results are the 
most prominent. In particular, we aim to answer the following questions: (1) 
Are the phase transitions of second order for the entire range of the anisotropy 
strength, as predicted by CVT, or is there a tricritical point separating lines of 
the second- and first-order transitions, as predicted by both MFT and EFT? 
(2) Is there a reentrant phenomenon in the second-order phase boundary, as 
suggested by the CVT results but not by MFT nor EFT? (3) Does the line 
of first-order transitions situated within the ordered ferrimagnetic region, ob- 
tained by CVT but not by MFT nor EFT, really exist? (4) Can the system 
display up to two compensation points, as predicted by both MFT and EFT 
(not investigated by CVT)? 
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2 Model and Monte Carlo simulations 



The model of the mixed spin-1 and spin-3/2 Ising system on the square 
lattice is described by the Hamiltonian 



h=-jj: sfsf - j:isfr - ^sfr, (i) 

(ij) i 

where = ±|, ±| for A ions, = ±1, for B ions, J < is the nearest- 
neighbor coupling parameter between the ions on A and B sublattices, and 
Da, Db are the single-ion anisotropies acting on the spin-3/2 and spin-1 ions, 
respectively. In this study we will consider the anisotropy of a uniform strength 
i.e., D = Da = Db. 

A simulated L x L square lattice consists of two interpenetrating sublat- 
tices, each one comprising sites. We consider linear lattice sizes rang- 
ing from L = 20 up to L = 200 with the periodic boundary conditions im- 
posed. Initial spin states are randomly assigned and the updating follows the 
Metropolis dynamics. The lattice structure and the short range of the interac- 
tions enable vectorization of the algorithm. Since the spins on one sublattice 
interact only with the spins on the other, each sublattice can be updated si- 
multaneously. Thus one sweep through the entire lattice involves just two sub- 
lattice updating steps. For thermal averaging, we typically consider N = 10^ 
MC sweeps in the standard and up to = 10^ MC sweeps in the histogram 
MC simulations [23l[2l], after discarding another 20% of these numbers for 
thermalization. To assess uncertainity, we perform 10 runs, using different 
random initial configurations. Then the errors of the calculated quantities are 
determined from the values obtained for those runs as twice of the standard 
deviations. 

We calculate the internal energy per site e = E/L"^ = {H) / L"^ and the 
sublattice magnetizations per site 




niB = 2{Mb)/L^ = 2{ 



B 



B 



(3) 



where (■ ■ ■ ) denotes thermal averages. The total magnetization per site is 
defined as 

m = {M)/L^ = {Ma + Mb)/L^ (4) 

Since for ferrimagnets m can vanish within the ordered phase at the compen- 
sation temperature, as an order parameter it is useful to define the staggered 
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magnetization per site as 

m, = {Ms)/L^ = {Ma - Mb)IL\ (5) 

Further, the following quantities which are functions of the parameters E 
or/and O (= M, Mg) are defined: the specific heat per site c 

the direct (O = M) and staggered (O = Mg) susceptibilities per site xo 

the logarithmic derivatives of (O) and (O^) with respect to (3 = l/ksT 

fl.o = ^ln{0^) = ^-{B), (9) 
the Binder parameter U 

and the fourth-order energy cumulant V 

V = l-^^. (11) 

3(^2)2 ^^^^ 

The above quantities are useful for localization of a phase transition as well 
as for determination of its nature. For example, the transition temperatures 
can be estimated from the locations of the peaks of the response functions c 
and Xo or, with a higher precision, from the intersection of the Binder pa- 
rameter U curves for different lattice sizes L. The first-order character of the 
transition can be judged from the presence of discontinuities in the magneti- 
zations and energy, as well as from hysteresis loops. Furthermore, the order of 
the transition can be reliably established by a finite-size scaling (FSS) anal- 
ysis. For instance, the energy cumulant V exhibits a minimum near critical 
temperature Tc, which achieves the value ^* = | in the limit L — > oo for 
a second-order transition, while l^* < | is expected for a first-order transi- 
tion [23l[2l] . Temperature-dependences of a variety of thermodynamic quanti- 
ties display extrema at the L-dependent transition temperatures, which at a 
second-order transition are known to scale with a lattice size as, for example: 

XO,ma.{L) (X L-^O'^O ^ (12) 

Dio,ma.{L) oc L'/'^o , (13) 
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D20,ma.{L) ^ L^'-o ^ (14) 

where vq and 70 represent the critical exponents of the correlation length 
and susceptibility, respectively. More precise locations of the extrema used in 
FSS can be obtained from histogram MC simulations [231121], performed at 
criticality for each lattice size. In the case of a first-order transition, typically 
a bimodality appears in the P{E) and P{0) histograms as L increases and 
the quantities (16|) — ( ITT]) display a volume-dependent scaling, oc L^. 



3 Results and discussion 

[Fig. 1 about here.] 

The critical temperature as a function of the anisotropy strength shown 
in Fig. [H was determined from the positions of the staggered susceptibility XMs 
peaks obtained by standard MC simulations for a fixed lattice size of L = 40. 
We checked that the peaks positions (the temperatures at which they occur) 
do not significantly change when larger L is used. For the case of Da = Db = 
we estimated the critical temperature with relatively high precision from the 
Binder parameter analysis, using the lattice sizes from L = 40 up to 200 
and = 5 X 10^ MC sweeps, as keTJlJ] = 2.362 ± 0.003. In line with 
the MFT P], EFT [1MB] and CVT [IS] results, the critical temperature 
decreases as the anisotropy is decreased from positive to negative values and 
eventually vanishes at the exact value oi D/\J\ = —2. However, in contrast 
to the CVT results, no signs of the reentrant phenomenon was observed near 
D/\J\ = —2. Namely, as temperature is increased from low values there is a 
single order- disorder phase transition for each value of D/\J\ > —2 and no 
phase transition below D/\J\ = —2. 

[Fig. 2 about here.] 

Furthermore, no indications of the additional first-order transition line 
within the ordered phase at low temperatures, found in the CVT study, were 
observed. In Fig. [2] we show anisotropy dependences of the staggered magneti- 
zation Tils and energy e, along with the respective response functions, XMs and 
c, shown in insets, at kBT/\J\ = 0.1. It is well known that discontinuities in 
physical quantities associated with a first-order transition tend to smear out 
in MC simulations for smaller lattice sizes. Therefore, in order to detect first- 
order transition features, we chose a relatively large size of L = 120 and also 
checked for possible loops in the curves obtained by increasing and decreas- 
ing D. According to [19], we should observe two anomalies associated with 
two phase transitions: the first-order one between two ordered phases and the 
second-order one between the ferrimagnetic and paramagnetic phases. How- 
ever, as evidenced in Fig. [21 only one anomaly can be observed and, therefore. 
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we conclude that there is only one phase transition which is from the ordered 
ferrimagnetic to the paramagnetic phase. This scenario agrees with the MFT 
and EFT predictions. However, the latter insist that in this low-temperature 
region, i.e., below the tricritical point J|, /c^Tj/l J|) = (—1.9730,1.1606) 
from MFT and {Dt/\JlkBTt/\J\) = (-1.9981,0.6175) from EFT, the tran- 
sition is first-order, which does not appear so in our simulations, since no 
discontinuities nor loops in the staggered magnetization and internal energy 
curves can be seen. 

Nevertheless, the respective spike-like response functions, which are some- 
what reminiscent of a first-order transition, prompted us to further explore the 
order of the transition by running more extensive histogram MC simulations 
for even larger sizes L and performing FSS analysis. We considered the lattice 
sizes up to L = 200 and the number of MC sweeps N = 10'', however, we did 
not detect any signs of a bimodal distribution in the energy or staggered mag- 
netization histograms at criticality that would signal a first-order transition. 
On the contrary, FSS of the quantities ([7]) — ([9]), shown in Fig. [3l indicates that 
the calculated susceptibility and correlation length critical exponents, 'JMs and 
uMsj respectively, agree well with the standard 2D Ising ones, i.e., 7/ = 7/4 
and I// = 1. In the susceptibility exponent dependence we can see that a linear 
regime in the log-log behavior is established only at L ^ 80 and, therefore, 
sufficiently large lattice sizes should be used. The second-order nature of the 
transition is also confirmed by the scaling of the minima of the fourth-order 
energy cumulant V (Fig. H]), which tend to the value V* = 2/3 for L ^ 00 
and apparently do not scale with volume, as expected for the second-order 
transition. Such thorough investigations were not done below fc^T/l J| = 0.1 
and we cannot rule out the possibility of the tricritical point existence at still 
lower temperatures but we deem it unlikely and presume that the transition 
remains second-order at all temperatures. 

[Fig. 3 about here.] 

[Fig. 4 about here.] 

Finally, we verified the MFT [15] and EFT [TS] predictions about the 
multicompensation behavior of the system. Namely, these theories predict 
the existence of two compensation points at the anisotropy values close to 
D/\J\ = —1.96. Such a behavior was indeed confirmed in our MC simulations. 
In Fig. [5] we plot the total magnetization vs temperature curves for several 
values of the anisotropy strength d = D/\J\. While for di = —1.95 there 
is no compensation point, for the values of ^2 = —1.958 and d^ = —1.97, 
there are two compensation points below the respective critical temperatures 
tc2{= kBTc2/\J\) and tcs- For d^ = —1.98, there is only one compensation point 
below tcA- The entire curve of the compensation temperature kBTk/\J\ as a 
function of the anisotropy strength D/\J\ is shown in the inset of Fig. [TJ 
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[Fig. 5 about here.] 



4 Conclusions 

We have studied the critical and compensation properties of the mixed 
spin-1 and spin-3/2 Ising ferrimagnet with a uniform single-ion anisotropy on 
a square lattice by Monte Carlo simulations. We focused on several contra- 
dictory results previously obtained by different approximative approaches. In 
particular, we checked if the system for sufficiently large anisotropy strength 
can display a reentrant phenomenon and additional first-order phase transi- 
tion within the ordered ferrimagnetic phase, as predicted by CVT. Neither 
of these predictions was confirmed. Instead, we found just one order-disorder 
phase boundary as a single- valued function of the anisotropy within the entire 
range of values, in agreement with the MFT and EFT results. Thus, while 
for the system on a cubic lattice the qualitative features predicted by CVT 
were confirmed by the MC simulations |22], they were not confirmed on a 
square lattice. On the other hand, our findings neither support the MFT and 
EFT predictions about the existence of a tricritical point at which the transi- 
tion would change to the first order one. Nevertheless, the multicompensation 
behavior with two compensation points observed within the MFT and EFT 
studies was reliably verified in the current MC simulations. 
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List of Figures 



The critical temperature /csTc/|J| as a function of tlie 
single-ion anisotropy strength D/\J\. The inset shows the 
enlarged picture in a low-temperature region with both the 
critical (fc^Tc/l J|) and compensation (fc^Tfe/l J|) temperature 
curves. [lO 

(a) The internal energy e and specific heat c (inset) per 
site and (b) the staggered magnetization and staggered 
susceptibility XMs (inset) per site, as functions of the 
anisotropy strength D/\J\, at kBT/\J\ = 0.1. [ll 

Scaling behavior of the maxima of the quantities ([7]) — OH]) 

for the parameter O = Ms i.e., XMs,max, DiMs,max 

and D2Ms,max, in a log-log plot at the critical point 

{Dj\J\,kBTj\J\) = (-1.9975,0.1). The slopes represent the 

values of the standard 2D Ising exponents ratios '^i/vi = 7/4 

and = 1. [l2 

The minimum of the fourth-order energy cumulant as a 
function of L^^. The solid horizontal line denotes V* = 2/3. |l3 

The total magnetization vs temperature curves for several 
values of the anisotropy strength d = D/\J\. There is no 
compensation point for di = —1.95, two compensation points 
for d2 = —1.958 and d^ = —1.97 and one compensation point 
for = —1.98. tci = ksTci/lJl, i = 1,2,3,4, denote the critical 
temperatures corresponding to the respective values of d^. |l4 
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Fig. 1. The critical temperature fesTc/] J| as a function of the single-ion anisotropy 
strength D/\J\. The inset shows the enlarged picture in a low-temperature region 
with both the critical (fcsTc/l J|) and compensation (/csTfc/l J|) temperature curves. 
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Fig. 2. (a) The internal energy e and specific heat c (inset) per site and (b) the 
staggered magnetization and staggered susceptibihty (inset) per site, as 
functions of the anisotropy strength D/\J\, at kBT/\J\ = 0.1. 
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Fig. 3. Scaling behavior of the maxima of the quantities ([7]) — ([9]) for the parameter 
O = Ms i.e., XMs^max, Dim,, max and D2Af,, max, in a log-log plot at the critical point 
{Dc/\J\,kBTc/\J\) = (—1.9975,0.1). The slopes represent the values of the standard 
2D Ising exponents ratios 'ji/uj = 7/4 and l/vj = 1. 
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Fig. 4. The minimum of the fourth-order energy cumulant 1^ as a function of L ^. 
The sohd horizontal hne denotes F* = 2/3. 
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Fig. 5. The total magnetization vs temperature curves for several values of the 
anisotropy strength d = D/\J\. There is no compensation point for di = —1.95, 
two compensation points for d2 = —1.958 and d^ = —1.97 and one compensation 
point for ^4 = —1.98. td = kBTci/\J\, i = 1,2,3,4, denote the critical temperatures 
corresponding to the respective values of di. 



14 



